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Abstract: We discuss fermion coupling in the framework of spinfoam quantum gravity. We analyze the gravity- 
fermion spinfoam model and its fermion correlation functions. We show that there is a spinfoam analog of PCT 
symmetry for the fermion fields on spinfoam model, where a PCT theorem is proved for spinfoam fermion correlation 
, functions. We compute the determinant of the Dirac operator for the fermions, where two presentations of the Dirac 
determinant are given in terms of diagram expansions. We compute the fermion correlation functions and show that 
^ they can be given by Feynman diagrams on the spinfoams, where the Feynman propagators can be represented by a 
discretized path integral of a world-line action along the edges of the underlying 2-complex. 
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1 Introduction 

Loop Quantum Gravity (LQG) is an attempt to make a background independent, non-perturbative quantization of 
4-dimensional General Relativity (GR) - for reviews, see [1-3]. It is inspired by the classical formulation of GR as a 
dynamical theory of connections. Starting from this formulation, the kinematics of LQG is well-studied and results in 
a successful kinematical framework (see the corresponding chapters in the books [1]). The framework of the dynamics 
in LQG is still largely open. There are two main approaches to the dynamics of LQG, they are (1) the Operator 
formalism of LQG, which follows the spirit of Dirac quantization or reduced phase space quantization of constrained 
dynamical system, and performs a canonical quantization of GR [4]; (2) the Path integral formulation of LQG, which 
is currently understood in terms of the spin- foam formulation [3, 5-8]. The relation between these two approaches is 
well- understood in the case of 3-dimensional gravity [9] , while for 4-dimensional gravity, the situation is much more 
complicated and there are some recent attempts [10] for relating these two approaches. 

A serious shortcoming of LQG and the spin-foam models has long been the difficulty of coupling matter quantum 
field theory (see the first two references in [3]), especially the coupling with fermions. It is still not clear so far about 
what is the behavior of the matter quantum fields on the quantum background described by LQG, and what are 
the quantum gravity corrections for matter quantum field theory. There was early pioneer works on coupling matter 
quantum field theory in canonical LQG [11] and in the context of spin- foam models in 3-dimensions and 4-dimensions 
e.g. [12-14]. And there was recent progress in [15], where we define a very simple form of fermion and Yang-Mills 
couplings in the framework of a 4-dimensional spin-foam formulation. Because of the simplicity of the fcrmion-coupling, 
it is possible for us to further analyze the detailed properties of the quantum fermion fields coupling to spin-foam 
quantum gravity. 

In the present article, we mainly discuss the fermion coupling in the framework of 4-dimensional Lorentzian EPRL 
model [6] ^ . The EPRL model in LQG is mostly inspired by the 4-dimensional Plebanski formulation of GR (Plebanski- 
Holst formulation by including the Barbero-Immirzi parameter /?), which is a BF theory constrained by the condition 

^Tho fermion coupling analyzed in this artiele can also be translated into Euclidean signature and implemented in the Euclidean 
EPRL-FK model and the model defined in [8]. 
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that the B field should be "simple" i.e. there is a tetrad field such that B — *{e A e). In the EPRL model, the 
implementation of simplicity constraint is understood in the sense of [16]. More importantly, the semiclassical limit 
of EPRL spin-foam model is shown to be well-behaved in the sense of [17, 18]. 

Our analysis of fermion-coupling in this work follows the definition in [15]. In Section 2, we review the regularization 
procedure of the fermion action on a 2-complex K. and discuss its formal continuum limit. We also show that there 
is a way to express the Dirac fermion action (more precisely, the Dirac operator) in terms of spin-foam variables, so 
that the Dirac action is coupled into the spin-foam amplitudes. In this way we describe the dynamics of the fermion 
quantum field theory on a quantum background geometry, which is described by spin-foam model. Moreover we define 
and discuss the fermion correlation function on spin-foams. In [12] it was mentioned that the (non-gauge-invariant) 
fermion correlation function vanishes on 3-dimensional spin-foam quantum gravity (the same thing also happens in 
lattice gauge theory, see the first reference of [12] and the references therein), so one should make a certain gauge- fixing 
in order to define the correlation function properly. In 4-dimensions we would find the similar vanishing result for 
(non-gauge-invariant) fermion correlation function, if there was local SL(2,C) gauge invariance. However in Lorentzian 
spin-foam model for pure gravity, a gauge-fixing has been implemented in order to make the vertex amplitude finite 
[19], such a gauge-fixing breaks the local SL(2, C) gauge invariance, and makes the fermion correlation functions well- 
defined. After that we discuss the PCT symmetry of the spin-foam fcrmions. The invariance under the inversions 
of charge, parity, and time simultaneously is believed to be a fundamental symmetry of nature. Here we define a 
transformation O for the gravity-fermion spin-foam partition function and its fermion correlation function, which can 
be viewed as a spin-foam analog of the PCT transformation in standard quantum field theory. Then a spin-foam PCT 
theorem is proved for the spin-foam fermion correlation functions, which states that the complex-conjugated fermion 
correlation function on a spin- foam background equals the correlation function of charge-conjugated fermions on a 
time-and-space reversed spin-foam background. This result is considered as a spin-foam analog of the celebrated PCT 
theorems proved for the quantum field theory on Minkowski spacetime [20] and on curved spacetime [21]. 

In Section 3, we continue the computation for the gravity-fermion spin-foam model. If the integrations of the 
fermionic variables are carried out, it results in a determinant of the Dirac operator on the spin-foam model. This 
Dirac determinant contains the information about the interaction between the fermion field and gravitational field. 
So in Section 3 we provide two representations for computing the spin-foam Dirac determinant in terms of diagrams. 

In Section 4, we compute the n-point correlation functions of spin-foam fermions. It turns out that the resulting 
spin-foam fermion correlation functions can be understood as coupling free-fermion Feynman-diagrams into the spin- 
foam amplitude, while each amplitude from a Feynman diagram depends on the spin-foam background geometry, 
which is summed over in the spin-foam amplitude. Because here we only consider the interaction between fermions 
and gravity, the Feynman diagrams coupled with the spin- foams are free fermion Feynman diagrams, which are 
completely factorized into Feynman propagators (matrix elements of inverse Dirac operator). And it turns out that 
the Feynman propagators can be expressed as a discretized path integration of a certain world-line action, where 
the world-lines are along the edges of the 2-complex underlying the spin-foam amplitudes. Our results confirm to 
some extend the early idea in [13] which proposes the inclusion of matter quantum fields by coupling their Feynman 
diagrams into spin-foam model. 

2 Gravity- Fermion Spin- foam Model and PCT Symmetry 
2.1 Definition of Spin-foam Fermion 

Given a 2-complex /C embedded in the 4-manifold M, we consider the discretization on the complex IC the classical 
Dirac action: 




(2.1) 
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where ^^{x) — 7^6^(0:) is the spinorial tetrad, e{x) — det(efj), and is the covariant derivative for Dirac spinor, i.e 



/I/ 
1 



D^4>^d^4> + i^K'Sij^ (2.2) 

here S^^ := 3:[7^,7''] is the Lie algebra generator of Lorentz group. 
We first consider the first term in the action, it can be written as 



5i := J / i^l'Di^ he' he"" A e^^euKL (2.3) 

which motivate us to make the following anzatz for the naive discretization of Si. We assign a fermion "^p^ to each 
vertex v E V{IC) of the complex JC, and make the following formally discretized Si [15, 22] 

5i ~ 2i ^ Ve V'h(e) 7^ ni{e) [G'eV'/(e) - V'Ke)] (2-4) 

ee_E(/C) 

where 6(e) and /(e) are respectively the begin and final point of e, n (e) = \ei^[v)e' (v)\ ~ ^(^)) ^ ^^i^ vector at 
the begin point of e, Ve is a 3-volume associated to the edge e, which can be viewed as the volume of the tetrahedron 
(polyhedron) dual to e ^, and 

Ge:=Ve^^''^'''^'' (2.5) 

is a SL(2,C) group element represented on the Dirac spinors. 

Let us check the formal continuum limit of Eq.(2.4), we consider a region f2, such that is much larger than the 
scale of a elementary cell (e.g. a 4-simplex), but smaller than the scale over which the fermion field and gravitational 
field (and their derivatives) change significantly. In this region f], along each edge, (ef — {d/dsiY) 

Ge,V/(e,) ^ [1 + As, e>t ]^Al^Su\ [V'fc(e,) + ef (9^^/')6(e,) 

~ + As, ef (i?^V')f,(e,) (2.6) 

Then at the region il, the formal continuum limit of Eq.(2.4) given by 

2i Ve V'fc(e) n/(e)7-^ As^ {D^^)b(e) = 2i ^ K V'&(e) ni{e)"f^ Ase e'^ e;^e;j(i:'„V)6(e) 

eCf2 eCO 

= '-Yo\{n)'^^"D^^P (2.7) 

where we have used the averaging formula: 

Vefijie) As, e^'e;^ = ^Vol(f^)5/ (2.8) 

To illurstrate this averaging formula: firstly because there is a very large number of edges (with all the possible 
directions) in the region Q, the left hand side of this formula is invariant under 4-dimensional rotation, thus is 
proportional to 6j . Secondly, if we take the trace of the left hand side, it gives the volume of the region il. 
Now we consider the discretized Si in Eq.2.4. The local SL(2,C) gauge transformations Uy act as follows 

Ge ^ C/b(e)GeC/7(^) V« ^ C/.V-. ^ ^vK^ nj ^ ^tr (jjUyj-^U-^) nj (2.9) 

The discretized Si is invariant under these gauge transformations. This can be seen by using the relation 

u.r'u-' = 7'^tr {^jUy-f'u-') (2.10) 



^Here we use r to denote both a tetrahedron and the center of a tetrahedron. We can make this notation because of their one-to-one 
correspondence . 



- 3 - 



Similarly, the complex conjugate term 5*2 in Eq.(2.1) can be discretized similarly 



- X! ['^e^/(e) -Aie)] Y nj{e) V'&(e) (2.11) 



2 

eeB(/C) 



while the mass term 5*3 is given by 



-mo J2 ^^-^vi^v (2.12) 



where ^Vi, is the 4-volume associate with v. 

As a result, the formally discretized action reads 

eeE{K) veV{K) 
= X! '5'e['(/'f,(e),''/'/(e),5e] + X! ^^i"^^] (2-13) 

eG-E(/C) t>ev(K;) 

where Ge is the representation of ge E SL(2, C) on Dirac spinors. 

So far the unit vector n^(e) are located at the begin point v of each edge e, and we assume n^(e) to be time- 
like and future-directed, which means that it can be transformed into (1,0,0,0) by a proper orthochronous Lorentz 
transformation. We make a parallel transportation of (e) from the begin point of e to a middle point r, such that 
n^(e) is transformed into n-^ {t) — Sq — (1,0,0,0), i.e. we consider a Lorentz transformation 

n^(e) Jtr {^jG,,^' G-,') = 5^ "/(e)^tr {^'G^rljG-,') = "/(e) Jtr {^jG-,WG.r) = 5'} (2.14) 

while 

A^j = ^tr (7^G-,Var) ^ (A-i)^^ = itr(7jG„.7'G-r) (2.15) 

niie)^^tT{^iG,rl'G-,')S^j (2.16) 

7^n,(e) = C7''G-,i55 - G„,7°G-,i (2.17) 

We define the Lorentz transformation Gyr as the spin-foam SL(2, C) holonomy g^r represented on the space of Dirac 
spinors. Recall that = (1, 0, 0, 0) is the unit vector orthogonal to all the face bivecors of the tetrahedron (polyhedron) 
r by the simplicity constraint [6, 16, 17]. So this definition states that the internal vector n^(e) coming from the tangent 
vector of the edge e is the normal of the tetrahedron (polyhedron) Te viewed in the frame at the vertex v. Therefore 
under this definition the fermion action is expressed by 



As a result 



then 



= -mo^V^i^.V. (2.18) 
We can also write Se in terms of Weyl spinors. In Weyl basis 

-(C) ^-(:;) 

A representation of SL(2, C) on Dirac spinors can thus be written as a tensor product between an SL(2, C) representa- 
tion on 2-spinors ^'^ and a SL(2, C) representation on dual complex-conjugated 2-spinors 9a' ^- Then the action can 

''Givon g G SL(2,C) a 2 X 2 complex matrix with a unit determinant, its representation on 2-spinors is given by h-> g^gS,^ , its 
representation on dual complex-conjugated 2-spinors is given by 6^/ i— > —g 8gi = —Sc' A'd'^jyi^^ ^ ^B' (from g^g,8^ ). 
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be written as 



= 2zK 



= -mo^V„ + em (2.20) 

where g^e is a SL(2,C) group element represented by 2 x 2 complex matrix with miit determinant. Here ^ denote 
complex conjugate and transpose^. In terms of spinor indices 

elfflff26 ^ <5s's(5i)^i'(52)^A(?i)-''(6)^ 

e\g\g^e2 ^ S^^'{gi)s^{g,)s'i'(e,)A{02)A' = 5BB'{gi)\{g2rk'{0i)^i02)^' ■ (2.21) 
Note that similar to Eq.(2.17), we also have 

CTi'A"/(e) '5?Cr|j,B(5^^_f,(e))^A'(5r,,6(e))^A and 0%^ = ^5b'B- (2.22) 

In the discretized fermion action, the 3- volume is a function of the spins jf and the intertwiner i^. More explicitly, 
we can express the vertex amplitude ^u[j/, ie, 5i;e]G-R in terms of Livine-Speziale coherent intertwiners = ||j/,f^e/) 
[23], whose labels (j/, fief) determines the geometry of a tetrahedron (polyhedron) [24]^. So we write the 3-volume Ve 
as a function of the labels (j/, n^f). For the 4- volume ^V^ of a 4-simplex, we can write it as a function of the boundary 
data (jy , fief) and SL(2, C) group element g^^ by the following procedure. Given a unit 4-vector (1, 0, 0, 0), we specify 
a so(3) sub-algebra inside the Lorentz Lie algebra (as a real Lie algebra), while we denote by the 4x4 matrices Li 
(i = 1, 2, 3) the so(3) generators in the Lorentz Lie algebra (in its vector representation). Given a tetrahedron dual to 
an edge e, its face bivectors are given by B^f — j fn\jLi. This bivector viewed from the frame of 4-simplex is defined 
by a parallel transportation B^f = QveBefQev = j fn-lfQveLiQev where g is the vector representation of g G SL(2,C). 
Given two triangles /, /' which don't belong to the same tetrahedron, we can write a 4-volume 

*Vy{fJ') =^ Byj,Byf, -f-^ jfjf,nljnl,f, < dveLtQev, Qve' LkQe'v >- (2.23) 

where the inner product -< -81,-62 ^uklBI'^ B^^ ■ The above way to write 4-volume depends on the choice of two 
triangles /, /', thus the expression of should be the above quantity ^V«(/, /') averaged by all the possible choices 
of /, /'. However in the large-j regime of EPRL model [17] or the model in [8], where the shape-matching condition 
is implemented, ^Vt, (/,/') doesn't depend on the choice of /, /', because the face bivectors can be expressed in 
terms of the wedge produces of cotriads [25] . 

We have expressed all the quantities in the fermion action in terms of spin-foam variables. Let's now couple the 
fermions to spin-foam quantum gravity. The spin-foam model for pure gravity on a 2-complex K. with a boundary 
graph 7 

ZGR(JC)f = / '^gveWAf[jf]GRT\A^[jf,ie,gve]GR ' Aj/.^e (2.24) 

where fj.jjj^ is a boundary state on the boundary graph 7 in its spin- network representation. In order to couple the 
Dirac fermion to gravity, we add an edge amplitude 

: 9ve:jf: *e] ' f^,jj.i^ 

where A}[j f] ^ Af[j f]GR, v4e[?/'fc(e), -0/(6), S^^e, j/, «e] = e*'^% A^[jg,ie, gve.tpv] ^ A,,[jg,ie, gve]GR e"^"' (2.25) 
The measure [DipyVip^] is defined by the Grassmann integral 

[D^.V^,] = n dV-^dC- (2.26) 

{v,a) 

■'Using spinor language, ^ for is defined by = ^a'Ai ^A' defined by {0'^)^ = 5'* '^Qa- 

^Tfien the sum over all the intertwiners in the spin-foam model will be a integral over all the unit 3-vectors n^j- satisfying the closure 
constraint [24]. 
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2.2 Fermion Correlation Functions 

With the gravity-fermion spin- foam model defined above, we consider a correlation function with a number of ipy and 
Tpy insertions: 



Spinfoam 



/ <igve TT Af[jf]cR TT Ay[jf, ie,gve]GR ' Aj/.^e ^ 

JSL(2,C) ^ 

X / [VipyVipy] ipy^--- V'i,„'0i'„+i ■ ■ • '0i)„+™ exp (^iSp [IC,jf, ie,gve,i^v] ) (2.27) 



The above definition of correlation function needs some explanations: Given that (1) the vertex amplitude Ay[jj:, ie, gvelcB. 
is invariant under the SL(2,C) gauge transformation g„e ^ ^vgve {^v €E SL(2,C)) at the vertex v; (2) the fermion 
action Sp [JC, j f,ie, gve,'4'v] is invariant under the SL(2, C) gauge transformation Gy^ ^ A^G^g, ipy i— A^ipy and 
ipy I— TpyA~^ {Ay and Gye are the SL(2,C) elements A„ and gye represented on the space of Dirac spinors); (3) 
the measures dgyg and dtpydtjjy is invariant under the SL(2,C) gauge transformation, we make a change of variables 
gve '-^ gye — ^vgve, i^v ^ i^v = Ayipy and il)y i-> ijjy — ^jjyAy^ of the integral (which doesn't change the integral at 
all), then use the gauge invariance of the vertex amplitude, fermion action, and the measures 



E 



Spinfoam 

- . . \ 1 

X 



[if]GR. TT Ay [jf , if, , gye]GR ■ f-r^jf 

JSU2.C) 



X J [V^^Vi^i] • • • ■ ■ ■ exp (iSp [lC,jf,ie, gt. ^v] 

= E / '^gve\\ Af[if]GB.\\Ay[jf,ie,gye]GR- 1-1,0 y- 

X j [D^yVi^J V^^^ • • • ■ ■ ■ "Pt+m exp {^Sp [JC, jf, ie,gve,^v] ) (2.28) 

Therefore naively one would have 

(^., ■ • • i^.^v^^^ ■ ■ ■ V'„„,„ )^ . ^ = (i^'y, ■ ■ ■ V^K,, • • • )^ . ^ 

\ / bpmioam \ / bpmioam 

= A,„„^, • • • A,„+„ (^y^ ■ ■ ■ i^yjy^^, ■ ■ ■ V'„„+,„ )^ ^ ■ ■ ■ K! (2-29) 

\ / bpinioam 

which results in that if wi, • • • ,Vn+m are different vertices and the correlation function is non-gauge- invariant, then it 
vanishes because the gauge transformations are independent at different vertices. The same result was mentioned 
in [12] in the context of 3-dimensional gravity, where it was suggested that one should either make a gauge-fixing to 
define the non-gauge-invariant correlation function, or instead consider the correlation functions of gauge invariant 
quantities, e.g. the fermion currents. 

However the above argument didn't take into account the regularization of the Lorentzian vertex amplitude [19]. 
If we consider pure gravity amplitude only, and our vertex amplitude Ay[j f ,ie,gyg\GR has the gauge invariant under 
gy,. I—)- Xygye, then the SL(2,C) integrals 



L 



TT dgye Ay[jf,ie,gye]GR (2.30) 

SL(2,C) e at V 

is divergent, because we can always choose an edge cq at v and a gauge At, = gyeo^ such that gye ^ 5Jeo5t'e (e 7^ eo) 
and gyef^ 1— >■ 1, then there is a redundant SL(2, C) integral Jd^yg^ which gives the divergence. The same argument 
for divergence also applies to the gravity-fermion spin-foam model by the gauge invariance of the fermion action 
Sp [)C, jf,ie, gve,^v]- The way to remove the divergence is firstly to choose an edge eo at each vertex v, and fix the 
gauge Xy = gy^^ at each vertex, such that gye H> .gtTeoSfe (e 7^ eo) and gyeg 1— > 1, so the integrand doesn't depend on 
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the variable g^^a for each vertex. Then the redundant integral J dg„eo should be removed at each vertex. It turns out 
that at least for the pure gravity spin-foam amplitude, the SL(2, C) integrals lead to finite result after regularization 
for a large class of spin- foam vertices [19]. 

Therefore the definition of fermion correlation function Eq.(2.27) should be understood in terms of the gauge- 
fixed/regularized spin-foam amplitude, where there is a edge eo at each vertex v such that the SL(2,C) group element 
gveo is gauge fixed to be the identity 1 e SL(2,C), and the integral of gyeg is removed at each vertex v. The previous 
argument for vanishing correlation function doesn't apply for the gauge- fixed/regularized spin- foam amplitude because 
the gauge-fixing breaks the SL(2,C) gauge invariance. Thus the correlation function Eq.(2.27) is not necessary 
vanishing. More explicitly we could write the correlation function in the following way instead of Eq.(2.27) 

{i'vi ■ ■ ■ V'i,„'0i.„+i ■ • • V'i.„+™)spi„foam 

<igveY[Af[jf]GB.Y[Ay[jf,ie,gve]GR- ff,jf,z^ 'X Y[ '5sL(2,C) (Sfeo ) X 

X y [PVi; 25-01,] V'l,! •• •V'«„'/'^;„+i •• + exp(iSF[IC,jf,ie,gve,i^v]^ (2.31) 

More detailed computation on the fermion correlation function will be given in Section 4. 
2.3 A PCT Theorem on Spin- foam 

We first briefly recall the notion of PCT symmetry for Dirac fermion on flat and curved spacetime. Given the 
fermion field operator ip{t, x) on Minkowski spacetime, its Parity-Inversion, Time-Reversal, and Charge-Conjugation 
are defined by 

P^{t,x)p-^ = j°ijj{t,-x) T^{t,x)T-^ = --i^-i^ip{~t,x) Cijj{t,x)C-^ = -ij^^*(t,x) (2.32) 

Thus the anti-unitary PCT transformation acts on the fermion field operator by (here we use ip* to denote complex 
conjugation and ■0*-^ to denote charge-conjugation) 

0'^(i,f) = {PCT)ij{t,x){PCT)-' = {~th'r{-t,-x) = /[^(„i^„2.)]t (2.33) 

Then the PCT theorem on Minkowski spacetime states that: Given the Minkowski vacuum state fl, a class of local 
fields ^{x) (can be composite fields) with m primed spinor indices, and an anti-unitary PCT operator {PCT) 

{PCT)^{x){PCT)-^ = <i>^{x) = (-l)™(-i)^$(-x)*, {PCT)n = n (2.34) 

where F = if $ is bosonic field and F = 1 if $ is fermionic field, the PCT theorem on Minkowski spacetime 
states that the complex conjugated correlation function of a number of fields equals the correlation function of the 
corresponding charge-conjugated fields [20], i.e. 

(f)|$f(xi)---$^(a;„)|0) = (^!|$l(a;l)•••$„(x„)|^^)^ (2.35) 

On Minkowski spacetime the map p : (t, x) h- !■ (— t, —x) defines an isometry preserving the spacetime orientation 
but reversing the time orientation. For a general (globally hyperbolic) spacetime Ai with metric gap and time 
orientation and spacetime orientation o = (T, SaP'ys)^ we let M be the spacetime with the identical manifold structure 
and metric structure as A4, but its time orientation and spacetime orientation is given by — o — (— T, Ea/s-ys)- Because 
there is no preferred vacuum state on a general curved spacetime, the PCT theorem on a general spacetime [21] is 
formulated in terms of Operator-Product-Expansion (OPE) coefficients: Given a local field $(a;) with n unprimed 
spinor indices and m primed spinor indices, we define its charge-conjugated field by 

$^(x) := i-in-if^{x)* (2.36) 
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Here x denotes the points on the spacetime manifold, and <i>'^(a;) is a local field with m unprimed spinor indices and n 
primed spinor indices. We suppose the a class of fields $^^(2;) (j = 1, • • • ,n) on the spacetime have the following 
OPE in short geodesic distance 

*ii,nyi)---<^(2/n)-E^>i..(2/i>--- ^M^^) (2-37) 

(3) 

as {yi, • • • , Un) approaching x, where the distributions c^j)^ xiv^^ ' ' ' > V") OPE structure coefficients. Then the PCX 
theorem on a general curved spacetime implies that on the spacetime Ai with the opposite time and space orientation, 
the OPE of the charge-conjugated fields ^j^*^^^ {j = ^, ■ ■ ■ ,n) is given by [2f] 

'J'S''(yi)--"fM''(y")-E^M,.(yir-- ,ynr <fg^(x) (2.38) 

U) 

whose OPE structure coefficients are the complex conjugation of c"J^ xiv^^ ' ' ' ' Vn)- 

In the formalism of spin-foam model, the quantization of spacetime M with metric g^p is formulated by a spin-foam 
amplitude Z{IC) on a 2-complex /C. The reversal of time orientation while keeping spacetime orientation unchange 
o I— )• — o can be formulated by simultaneously reversing all the internal edge orientations in the complex JC. The 
reason is the following: Given a spacetime {A4,gap) with space and time orientation o = {TjEap-ys)- AH the oriented 
orthonornial frames are given by e", satisfying 

ff„/3e?e^ = Tjij e^^c^T > e^^^s ef > (2.39) 

The oriented orthonormal frames form the frame bundle F{A4) over A4 whose structure group is the proper or- 
thochronous Lorentz group. However for the spacetime {AA,gap) with opposite space and time orientation — o = 
(— T, £^^75): its frame bundle F{M) is naturally isomorphic to F{M) by the map 

/ : F{M) F(M), ^ -e? (2.40) 

We consider the spin-foam model defined on a 2-complex /C dual to a simplicial complex imbedded in the spacetime 
manifold. Given a tetrahedron r in the simplicial complex, all the area bivectors of the tetrahedron are orthogonal 
to a unit internal vector = (1,0,0,0) by the simplicity constraint [6, 16, 17]. And this vector is given by the 
tangent vector of the edge e dual to the tetrahedron r up to a proper orthochronous Lorentz transformation, i.e. on 
the spacetime manifold Ai with orientation o 

aA[jn^ = e°'ei {a > 0) (2.41) 

where A^j is a proper orthochronous Lorentz transformation. However if the spin-foam model is build on the spacetime 
manifold A4 with orientation — o (in another words, if the spin-foam model is a quantization of the spacetime structure 
{A4, gap, —o)), the previous co-frame field changes into — e^, thus one has to change the previous tangent vector 
e" into — e", i.e. reverse the edge orientation, to keep unchange. Note that has to be fixed to be (1,0,0,0) for 
the definition of spin-foam model. 

We first formulate the PCT invariance of the spin-foam fermions in the following formal way: We define an 
anti- unitary PGT operator 6, such that it acts on fermion field operators by 

e^pj-' = i-ih'^'i^l ei>J-^ = (-*X7°7' (2.42) 

which are charge-conjugate field operators. If we consider the fermion action Sf{}C) (defined on a complex /C) is an 
composite operator from the fermion field operator ■0^, and ijjy, Sp is PCT invariant in the sense that 

0Sf{IC)9-^ = Sf{IC-^) (2.43) 
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It is indeed the case. We consider the bihnear form 2iVe'4>yGye'y^Gev''4'v' appearing in the discretized fermion action, 
by using {7^,7*^} = 0, (7°)^ = (7^)^ = 1 and (7°)^ = 7° = (7°)*, as well as the assumption is anti-unitary 

= -2zK(-i)V'*7°7'G:,7"G:„,(-z)7'7°V^!,' = -2iVei^WGl,fG:,,fi^l, = -2iV;V'^G*„7°G*,,^t' 

= 2zye^„'G„.e7°Ge„V. (2-44) 

where \vc have treated 14 and G^e as c- numbers (gravity as external field), and we also used the fact that 7^ commutes 
with G and the relation 7°G*7'^ = (G*^^)^ = (G~^)*. Here we have shown that the action of 9 on the bilinear form 
interchanges the vertices v and v' . Hence under this transformations of variables, S'e['!/'6(e)) V'/(e)] transforms to 

SebPb{e),tpf{e)] ^ '^iVe [V'/(e) Gf(^e)Tel°Greb(e)1pb{e)-'tPb{e)Gb{e)T,7°Gr,f(e) tpf{e) 

= S'e-i[V'6(e-i),V'/(e-i)] (2-45) 

where we use that 6(e) = /(e^^) and /(e) = b{e^^). On the other hand it is easy to see that the mass term 
Sy = —mo^Vvtpv''Pv is invariant under 6. Hence we obtain the PCT invariance of the spin-foam fermion in the sense of 

6SF{lC)e-^ = Sf{1C-^). (2.46) 

Now we consider the gravity-fermion spin-foam model. We define the following transformation O of the spin-foam 
amplitude, which is an analog of PCT transformation: 

Definition 2.1. Given a 2-complex K with a boundary graph 7, and the gravity-fermion spin-foam 

\jf\GRYYAy\jf,ie,gye]GR ■ fi,jf,ie X 

./,ie-^SL(2,C) f 

j[Vi;yV'^^] exp(^iSF[IC,jf,ie,gve,i^v]) (2.47) 



X 



where f-y,jj,i^ is a boundary state in its spin-network representation, we define a spin-foam analog of PCT transfor- 
mation O by 

• 9 reverses the orientations of all the internal edges in the complex K,, i.e. 9 : /C i-)- 1C~^; 

• 9 changes all the gravity vertex amplitudes ^t,[i/, ie, fftielci? into their complex conjugates Ay[jf,ig,gye]Gj^; 



The boundary state f'y,jf,i^ transforms into its complex conjugate f^^jj^i^; 



• 9 changes the fermion action on the exponential iSpiK,) into 

e [iSpilC)] = -iSpilC-^y, (2.48) 

• For each fermion ipy at a vertex, Q : ipy t-^ ip^ := {—i)j^4'l = {~i)7^1^''Pv (charge-conjugation), or in terms of 
Weyl spinors: 

e: 

The second and third transformations are motivated by the anti-unitarity of the PCT operator. Wc consider 
heuristically the pure gravity spin-foam amplitude as a physical inner product between some certain in-state and 
out-state fin and /out, which describe the boundary data of quantum gravity 

ZGR{lC)f = (/«„t , fin)pHys (2-50) 
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Then licuristically the PCT transformation reverses the in-state and out-state by 



^ • {fout ■) fin)pfiys ^ {Pfout 1 ^fin)pfiys {fout •> liti) Phys (/'in i foutjpfiyg 



(2.51) 



where we implicitly use the fact that the path integral measure Ag^e is real, since the SL(2, C) Haar measure dg can 
be written explicitly by 



d5 = 



d^d/3*d7d7*d(5d(5* 
W 



9 



a (3 
7 6 



(2.52) 



which is manifestly invariant under g g* . 

Here we argue that the PCT transformation of pure gravity spin-foam amplitude 



ezGRiic)f ^ ZGRiic-': 



f 



(2.53) 



is a spin-foam analog of the spacetime A4 with reversed time and space orientation — o = (— T, £^^^5). We consider 
the semiclassical behavior of the pure gravity spin-foam amplitude. For a given 4-simplex, the corresponding vertex 
amplitude Ay[jf,ie,gve]GR can be represented in terms of Livine-Speziale (LS) coherent intertwiners [23], i.e. use LS 
coherent intertwiner \ \jf,nf) for the SU(2) intertwiner i^. Thus we can denote the vertex amplitude by 



^v[jf,nf,gye]GR 



(2.54) 



where ri/ e S*^ is a unit 3- vector associated to each triangle/face /. The large-j asymptotics of the vertex amplitude 
^v[jf,n'f,gye]GR was studies in [17]: If jf goes to be large 



(2.55) 



where N± are independent of j/, and i9/ is the extrinsic angle between the two tetrahedra sharing the triangle /, i.e. 
for two tetrahedra r, r' sharing a triangle /, n-' [t] and {t') are two unit 4-vectors respectively orthogonal to all the 
face bivectors of r and r'. The extrinsic angle d f is defined by cosh??/ = —rjijn'^{T)n'^{T'). Then 



^ijf'^Ajf^nf] = SRegge\jf,nf] 

f 



(2.56) 



is the Regge action of discrete GR with triangle areas Af — SirGh/Sjf {/3 is the Barbero-Immirzi parameter). 

The vertex amplitude can also be written in the holomorphic representation [18]. Given the boundary graph 7 of 
the vertex amplitude, and for a given link I, we associate a complexifier coherent state t/;* (hi) [26] where 



(2.57) 



where Ms a dimensionless classicality proportional to £p, and gi E SL(2, C) is the complexified phase space coordinate 
gi = hie'^'^^^^p^ . These coherent states (with all possible gi) form a over-complete basis of the Hilbert space L^(SU(2)). 
The representation of the vertex amplitude on the complexifier coherent states is called the holomorphic representation 
of the vertex amplitude, denoted by ^i, [ff/, (7t,e]Gi? and is regarded as an amplitude with boundary data 5/ (a link I 
of the boundary graph uniquely corresponds to a face /). In the holomorphic representation and as the limit £p 
[18, 26] 



Av[gf,gve]GR^'^ K [j/ , n/ , g^e] exp 



J/ 



-^^3 lis 
f 



exp 



E 



2 ^ 



(2.58) 



where Af is the area of the triangle / evaluated at the phase space point g/, and — f^'d is the extrinsic angle 
evaluated at the phase space point gj. From this we see that as the limit Ip — > 0, the large spins jf ^ dominate 
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the contributions of [ff/, (?t,e]G-R- We then consider the large spin contributions and insert in the large-j asymptotics 
of the vertex amphtude Eq.(2.55) 



p=±i jf 



exp 



Af_yt 



(2.59) 



where the factor exp[i jf {p/S-df — ^f)] is a rapid osciUating phase as jf is large, so that the sum over jf suppresses 
unless the boundary data coincide with fi-df or —P'df- Without losing generality, we assume the boundary data 
^ = Pdf, thus in Ay[gf,gvf.]GR the terms with p = +1 are preserved in the sum, corresponding to e*'^"=99=, while the 
terms with p = — 1 suppress in the sum over jf. Therefore 



Av[gf,gve\GR exp 



exp 



PR) 2^Z^\PbJ 2 



(2.60) 



If we make the PCT transformation on the spin-foam amplitude Ay[gf,gye]GR, the reversal of the internal edge 
orientations doesn't affect the amplitude while the complex conjugation gives 



0^u[5/,.9^e]Gfl = A^[9f,gye]*GR 



exp 



Of 



exp 



E 



Jf 



IIP 



E 



(2.61) 



where effectively the terms corresponding to e~''^^='"'= is preserved while the terms corresponding to e''^^"'"' are 
suppressed. The flip from e^'^^^sf" to e~'^^"-^^<''' corresponds to the flip of the extrinsic angles from i9/ to — t?/, which 
is a discrete analog of flipping the sign of the extrinsic curvature by flipping the timc-oricntation. 

When we consider the fermions in spin- foam model, in terms of the semiclassical limit, the fermion coupling 
giSF should couple with the e*'^''=s3'=-terms (with positive coupling constant) on the spacetime manifold A4. It means 
that one should specify the boundary data such that e~*'^''=9»'-terms are suppressed, in order to reproduce the right 
equation of motion of graviy-fermion system on M with a positive coupling constant k — SttG. However when we flip 
the time and space orientation, e^'^^^ss" -terms are suppressed while e~*"^'*=!'9'= -terms representing the spacetime(s) A4. 
Then the fermion coupling has to flip sign as e~*'^^ to couple with e~*'^^'=9s% in order to keep the coupling constant 
positive and a right equation of motion. This gives an argument^ for the minus sign in the result of Eq.(2.48). 

We consider a fermion correlation function of charge-conjugated fermion fields ip^ = Oip^O"^ and = 0ip.^j6~^ 
on a PCT transformed gravity-fermion spin- foam amplitude, i.e. ( Spinfoam(A^) denotes the spinfoam analog of the 
spacetime M ) 



Spinfoam(A^) 



E / '^9veY[Af[jf]GRYlAy[jf,ie,gyXGR- fj.jf.t^^ 

Jsu2.c) r , ' ■ 



First of all we note that the the fermion measure is invariant under charge-conjugate 
We then make a change of variables for the fermionic integration 



(2.62) 



(2.63) 



(2.64) 



^The semiclassical arguments in this paragraph and the previous paragraph are heuristic arguments supporting the definition of the 
PCT transformation 0. A more precise argument relies on the large-j asymptotic analysis for the gravity-fermion spin-foam model on a 
large number of 4-simplices, which will bo studied in the future publication. 
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which doesn't change the integration. Then 

\ ^ ^ I Spinloam(Aa) 

X Jpr,T^rv]rv,---rvy:^+,---rv„+^ exp(^-iSF[ic-\jf,ie,gve,^',]) (2.65) 

Again we consider the biUnear form 2iVeil>'yGve'y^Gev''4''y' appearing in the discretized fermion action, 

= 2zK(-«)^Ah"7'G.e7"Ge„'(-*)7'7¥l' = 2iKV'h°G.e7°Ge.'7°V^I' = 2iVeijlGiyGl,J^l, 

= -2z^eV^:,G:,,7°G:>: (2.66) 

where we used the fact that 7^ commutes with G and the relation j^Gj^ = (G~^)^, as well as the anticommutativity 

of V'j;- 'S'e[V'6(e)> V'/(e)>S'fe] ^^^^^ 

SeWb{e)^^'fie)'9ve] = -2iVe G}(^)^^7°G;^(,(^)V6(e) " V'6(e)<^6(e)re7°G;^/(e) V'/(e) 

= S'e-i[V'6(e-i),V'/(e-i),fi't;e]* (2.67) 



And again for the mass term 



Therefore we obtain that 



sM = -mo^v.v'.V'; = mr.r, = sm* (2.68) 

SF[lC~^,jf,ie,9ve,'ipv] = Sp [IC, j f ,ie, Qvc'^v]* (2.69) 



Then as a result 



/-fisL(2,c) y 



3f 

X 



= (^i-^i • • • V'^„V'^„+l • • • ^v^+Sspin{o.MM) (2-70) 

We summarize the result as a spin-foam analog of PCT theorem: 

Theorem 2.1. (PCT theorem for Spinfoam fermions) 

The correlation functions of fermions on the spin-foam analog of a spacetime A4 with a certain time and space 
orientation o = (T^Safj-ys) equals the complex conjugated correlation functions of charge-conjugated fermions on the 
spin-foam analog of the spacetime M with an opposite tim,e and space orientation —o = (— T, Eq,^-^^), i.e. 

{i^Z ■ ■ ■ ^vJZ,. ■ ■ ■ , . , = {^v. ■ ■ ■ i^vjv..^ ■ ■ ■ V'.„,^)* . ^ (2.71) 

\ ^ ^ I SpmfoamyM) \ ' Spznfoa'm(A4) 

3 Determinant of Dirac Operator on Spin-foam 

To simplify the formula, now we consider a massless chiral Weyl fermion on the complex /C 

5'e [^b(e) ' 0(e))5r, 6(e) '.9re /(e)] = 2^14 [cl(e)3l,fc(e)5'r,/(e)^/(e) - 4(e)ffI„/(e)3r,fc(e)C6(e) 

SF[^v,gve] = Y Se['lpb{e),'^fie), grMe)' 9T,f(e)] = ^i^lPv.v'^v (3.1) 
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where the Dirac operator has the foUowing matrix element 

and all the other matrix elements vanish, where we see that p is an anti-Hermitian matrix 

= -p^y (3.3) 

We define the Grassmann path integral measure by 

^mK,^^] := X{^{ivfmi)imvfml)2- (3.4) 

V 

Obviously the path integral of the fermion action gives the determinant of the Dirac operator (Dirac determinant) 

det^[g.e] = j e'^-K">«"=l = j i^MK,^^] e" ^i^--''^'' (3.5) 

3.1 Polymer representation 

We briefly recall the polymer representation of the Grassmann Gauss integral [27] . Given a matrix Q such that 

= MAj - K,j (3.6) 

where Mi are diagonal contributions and Kij are off-diagonal contributions, so the diagonal elements of Kij are 
assumed to vanish Ka = 0. In terms of the Grassmann Gauss integral: 

det Q = j drJidryi • • • dT^j^di^N ^' m^^.^^+y.^^^ 5j.K.,r,, ^3 

This determinant can be represented in the following way: we first draw N points as a lattice representing the indices 
I = 1, • • • ,N . Then we draw all possible polymer diagrams using the building blocks in Fig.l, such that every index 
point has precisely one incoming and outgoing line (a monomer is counted both as a incoming and outgoing line). We 
denote the set of all possible polymer diagrams by V and denote a polymer diagram by z. 

Monomer: C^^^^ ~ 



Polymer < — K- 

line: A 



1-3 



Figure 1. The building blocks of polymers: a monomer at a single index point and a polymer line connecting two different 
index points 



For each polymer diagram we write down its contribution for det Q 

^ (_l)Numbor of Polymer Loops "Q "Q j^^^ ^3 

i j,k 

where the index point i is attached by a monomer Mi and the index points j and k are connected by a polymer 
line Kjk, while the polymer lines forming closed polymer loops. Then the determinant detQ equals the sum over all 
possible polymer contributions: 



det Q = ^ 



(3.9) 



Now we consider the determinant of our Dirac operator p. The diagonal elements of P are zeros, thus the 
monomer contribution is not allowed. For the off-diagonal elements: 

P{b{e),A'y,{f{e),A) = ^Ve ^B' B^eb{e))\'{9eS(e))\ P{f{e),A');{b{e),A) = (^B'S (5e/(e))^' (ffebCe) )^ (3-10) 

We draw the lattice of all the index points j = {v, A), which is V{)C) x Z2. Then we draw all possible polymer loop 
diagrams using polymer lines Pij , such that every index point has precisely one incoming and outgoing line. For each 
polymer diagram z gV 

= (•_;^^Number of Polymer Loops 

^ (_l^Nuxnber of Polymer Loops "Q [-If {2V,)5s, j,{g, {gev)\ (3-11) 

(v' ,A')-(v.A) 

where = 1 if f ' = b{e) and v = /(e) otherwise e^iy = 0. Finally the Dirac determinant is represented as a sum 
over all polymer contributions: 

det^=^i?^ (3.12) 

For the massive Dirac fermion, the lattice of all the index points is V{K) x {1, 2, 3, 4}. And we not only need to 
consider the polymer lines, but also need to consider the monomers iMj-=(„_Q,) = imo^Vy. For each polymer diagram 
z gV we have 

= (-1)^™^^'- Loops -Q-j|^^-Q(_^^.^) 

^ (_l^Number of Polymer Loops (,^^4-^^^ ^_-^y,,^ ^^y^^ [G,,,J°G,,] (3.13) 

And 

det(^ + iM) = Y^ (3.14) 

3.2 e-loop representation 

Here we compute the Grassmann path integral in a more explicite manner. Because of the Grassmann variables, on 
each edge the exponentiated fermion action has a 9-term expansion 

e^S- = 1 - 4(g) A(e),/(e)0(e) + 1^ (cl(e) A(e),/(e)?/(e)) 

1 - 4(e)^/(e),6(e)^6(e) + ^ (4(e)^/(e),6(e)^6(e 

= ^l(e)Pb{e)J-{e)if{e) - C}(e):P/(e),6(e)Cb(e) + A(e),/(e)^/(e)) (C/(e)^/(e),b(e)^6(e 

+ 2 \^b{e)Pb{e)J{e)^f(e) ) +2 [^f{e)Pfie)Me)^b{e) ) 

"2 (4(e)^Ke),/(e)^/(e)) (4(e)^/(e),6(e)^b(e)) ~ ^ (4(e)^/(e),6(e)^6(e)) (cl(e) A(e),/(e)^/(e)) 

+ \ {^l{c)PHe),f{e}i.f{e)) ^,,-^P f (e) Me)^b{e)) (3.15) 

Given a Grassmann integration at a vertex v. it only affects the e''^='s with the e's connecting to v. We have 



/ 



d(c.)'d(4)id($,)2d(^t)2 n n e^^^' (3.16) 

e,b(e)—v e' ,f{e')=v 

We know that the only nonvanishing Grassmann integral is 

J d^de i^T'^^^"" ^ J dcM(ct)^de'd(ct)^ ^^'F'^^e = (3.17) 

Therefore there are a few examples for the possible contributions for the integral at each vertex v 
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Example 1. Consider two outgoing edges 61,62 and two incoming edges 63,64, v = b{ei) = 6(62) = /(es) = f{ei), 
we have a nonvanishing integral 



J dC.dCj 



(3.18) 



PvJ{ei)A'A^f(^ei) ^f(e2)Pf{e2),vB'B\ [^b{e3)Pb{e3),vC'cPvMe4)D'D^b{e4.) 
PvJ{ei)A'A^f(^ei) ^ f{e2)P f{e2),vB' b\ [Cb{e3)Pb{e3),vC'cPv,b{e4)D' D^bie^) 
= ~ (^ne2)^b(e3)Pf(e2),vB'BPb(e3),vC'C£'^'^^ (^^ ^ PvJ{e^)A>APvMei)D'D^,j{e^)£,b{ei)) 

Insert in the matrix element of Dirac operator 

Pb(e)J(e)A'A = 214 ^B' B{Seb(e))\'{9ef(e))\ Pf{e)Me)A'A = -214 SB'B{9ef{e))^A'{9eb{e))^A {3.19) 

the above integral equals 

= -IGVe^Ve^Vg^Ve^ ~^ f(e2)^b{e3)9 e2f (e2)^B'9 esb{e3)^C'^E' E^F' Fge2v%9e3V^C^^'^ 

„ \.A'D'- G' ^ H' c r „ G „ H (A gD ' 

^ 9eiv A'9e4V D'^G'GdH' Hgeif{ei) A9eib{e4,) DS/(ei)?fe(e4) 

Example 2. Consider two outgoing edges 61,62, h{e\) = 6(62) = v, there is a nonvanishing integral 

^v"^ Pv,f(ei)A'A^f(ei) ^v^ PvJ(ei)B'B(f(e,) ?/(ei)'^ Pf(ei),vC'C^v C/(e2)^ P f (62) ,v D' D^v 
Cv^ Iv^ PvJiei)A'APvJ{ei)B'B^f(^ei)^fiei) ^f{ei)^ 0(e2)^ Pf{ei),vC'cPf{e2),vD'D^v^v 



(3.20) 



^ Pv,fiei)A'APvJ{ei)B'B^f(e,)Cf{ei) '^/(ei)'' 0(62)"^ Pf{ei),vC'cPf{e2),vD'De' 



D' i 



-CD 



= -8Ktl42 



(5lit,5ei/(ei))A'A(5lii.ffei/(ei))B'BC/(ei)?/(ei) 



^^■{e^)i^{e2)i9l^f{e^)9e^v)c'c{gl^f(^^^)ge2v)D'De 



CD 



= -8K'K2 



A gB 



£AB^f(e)^ 



^'f{ei)^f{e2)(9l^f{ei)9eiv)c'ci9l^f(^e2)9e2v)D'D£ 



CD 



(3.21) 



Example 3. Consider an outgoing edges 61, and an incoming edge 62, &(6i) = /(62) = v, there is a nonvanishing 
integral 

j d^vd^l 2 ^v^ PvJ(ei)A'A£,f(ei) ^v^ PvJ{ei)B'BS.f(ei) ^f(ei)^' P f (ei) ,vC' dv ^6(62)^ P f (€2) ,v D' D^v 



^ Pv.f{ei)A>APv,f{ei)B'B^f(ei)^f{ei)\ [^/(ei)*" ^/(e2)^ Pf{ei),vC'cPb(e2),vD'DS 
-^V^iVe2 {9liv9eifiei))A'A{9Uv9eif{ei))B'B^f(ei)^f{ei) 
X(-1) ^nei)^f{e2)(9l,fl^e,)9eiv)c'c{9l2bie2)9e2v)D'De'^'^ 



-CD 



-CD 



[^'fiei)^ne2)(9l^f{e^)9eiv)c'c{9l^b{e2)3e2v)D'D£' 

Example 4. For each single outgoing edge e, v = b{e), we have the integral 

J d^vd^l ^ ?t,^ PvJ{e)A'ACf(e) PvJ{e)B'BCf(^e) ^f(e)^ Pf(e),vC'cCv ^/(e)^ P f {e) ,v D' d(,v 

= -^^e [£^B?/(e)^/(e)J [?/(e)?/(e)^C"D' 

Example 5. For each single incoming edge 6, u = f{e), in the same way 

J d^vd^l ^ PvMe)A'ACf{e) ^v^ Pv,b{e)B' BS,b{e) ^b{e)'^ Pb(e),vC'civ ^b{e)^ Pb{e),vD'D^v 



(3.22) 



(3.23) 



£AB^b{e)^b(e) ^b(e)^b{e)^C' D' 



(3.24) 
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All the contributions of the determinant dot p can be obtained by the integrals for each vertex, similar to the 
previous examples. And they can be represented graphically: 

• We draw an arrow-line for each ^lPv.v'£,v in Se and associate to its source and ^ to its target. The arrow-line 
(Fig. 2) represents Pv,v'A'A with A' at its source and A at its target (however in the follow graphic representation 
we often ignore the A', A label for the arrow-line, in order to simplify the graph). 

A' *A 

Figure 2. A fermion arrow-line. 



• For each edge we represent by Fig. 3, where we have a minus sign (—1) when orientation of a fermion arrow 
line in Se coincides with the orientation of the corresponding edge e. The parallel double-arrows correspond to 
the terms {^^P^)^ in the expression of Se Eq.(3.15). 



^s^= (i+^ + =:)(i+^ + t=) 




(-1) ► 

< < ' < 



Figure 3. All possible contributions from an edge, see Eq.(3.15) for the correspondence. 

• We assign the weights to the arrows and double-arrows (Fig. 4) 

Figure 4. The weights for arrow and double-arrow. 

• We represent e^^ (resp. e'^ ^ ) by a contractors connecting two targets (resp. two sources) of two arrow-lines 
(Fig. 5). Note that the e-contractors not only can contract the arrow-lines from different edges, but also can 
contract the double-arrows from a single edge. 

• Each vertex v £ V{)C) must choose precisely two incoming and two outgoing fermion arrow-lines from the edges 
connecting v. These fermion arrow-lines are contracted by two e-contractors. There are 7 types of contractions 
shown in Fig. 6. Note that it requires precisely the same number of incoming and outgoing arrows because one 
need precisely two ^ and two ^ to make the integral nonvanish. 

Type 1: Four fermion arrow-lines are from four different edges; 
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Figure 5. A e-contractor connects two sources or two targets. 




Typel. I I Type 2. 



Type 3. ^ Type 4. 



■ € ■ 





Type 5. 



Type 6. * ! !* Type 7. _ * 

Figure 6. The typical contractions for eacli vertex. 

Type 2: Two arrow-lines are from two different edges, but a double arrow is from another single edge; 
Type 3: A pair of opposite oriented arrows from a single edge, and two arrows from two different edges; 
Type 4: Two pairs of opposite oriented arrows from two different edges; 

Type 5: A double-arrow and a single arrow from a single edge, and another arrow from another edge; 
Type 6: two double-arrows are from two different edges; 
Type 7: two double-arrows are from the same edge. 

The contraction at each vertex makes the arrow-lines form close loops. There are two type of loops, i.e. nontrivial 
loops and trivial loops, see Fig. 8. We call these loops the "e-loops" because the neighboring edges have opposite 
directions. Note that a double-arrow from a single edge only can form a trivial e-loop by the above contraction 
rule. And each £-loop must contain even number of arrow-lines by construction. 
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I* 1 

I 



Nontrivial Loop Trivial Loop 

Figure 7. A 4-gon nontrivial e-loop and a trivial e-loop. 



• We denote a n-gon nontrivial e-loop by i„ {n > 2 is even) and a trivial e-loop by T. Each trivial e-loop 
contributes 



T = -2VhABe 



AB 



(3.25) 



Let's consider a non-trivial e-loop, for example a 4-gon e-loop with vertices Vi^V2,v^,V4^ (cyclic ordered) and Vi 
is a source node for two arrows 



u = j deide2de^d^4 (ef Pi2A'a^^) (ef p^2B'b^!^) Pzac-c^x ) (ef Piad'dS.^ 

= e""'^' Pi2A'A{-e'^'')P:i2B'Be^''^' p:iiC'ce^°PuD'D 
- (-1)16T4,K.K3K,(-1)^^^ xe^'^' [^£'£(ge,.J A'(.9e,.J^A 



.AB 



SF'Fi9e2V3)^B'{9e2V2Y 



.B'C 



= (-l)16l^eiV;,K3K4(-l)'^*tr [{9e,v^£9lrJi9erV2e9l2V2)i9e2V3^9lsvs)i9e3V,S9ly,)] (3.26) 

where £l„ is the number of arrows whose orientations coincide to the orientations of the associated edges. 
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Figure 8. A 4-gon nontrivial £-loop. 



Generalize it to n-gon e-loop L„ (choose vi to be a source node, Vi , are cyclic ordered) 

n 

i=l 

where f/e.e' is the SL(2,C) holonomy from the middle point to the middle point Tg/. 

• We draw all possible close e-loop diagrams {Tp} by using the possible edge contributions (Fig. 3) and the 
possible vertex contributions (Fig. 6). For each e-loop diagram, it consists a certain number of trivial e- loops 
T(i), i = 1, • • • , |7rp|, and a certain number of nontrivial e-loops L„(i), i = 1, • • • , I^CffIi where Tr^ and C^p 
are respectively the sets of all trivial and nontrivial loops in F^. Therefore we can write the determinant of the 



lV2^9e2V2} 



{9e 



_i£5e 



Ji9e„-,vM„vJ (3-27) 
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Dirac operator is a sum over all possible loop diagrams 

ITrpI |£r^| 

det P = J2 Yl Til) Y[ Lnii) (3.28) 

{r_p} i=l i=l 
^ ^ ^ 2^Numbor of Loops ^Number of Trivial Loops ^ ^ 

{r^} {^} 

Y\ (-l)^^"tr ige„v,£fe^v,){9e,V2£gl^v2) ' ' ' (5e„_2«„-i e5e„-i«„-i ) (5e„_i«„e5e„«„ ) 
{Nontrivial Loops} 

where {—J'} denote the set of all the arrows in the e-loop diagram Tp, so means the product of all the 

Ve associated with the edges corresponding to all the arrows, and note that for an edge there can be two arrows. 

One could follow the following two steps to construct each term in the sum Eq.(3.28) 

Step 1. We first ignore all the trivial £- loops. All nontrivial e-loops can be constructed by using the first 4 terms of 
Se in Fig. 3 and the type 1-5 vertex contributions in Fig. 6 (ignoring trivial loops). But one should make sure that 
each vertex has 4 fermion arrows (2 incoming and 2 outgoing arrows) or has 2 fermion arrows (both incoming 
or both outgoing). 



Figure 9. A simple diagram with 2 nontrivial loops. 

Step 2. We add trivial e-loops such that all the vertex has precisely 4 fermion arrows, two of which are incoming and 
two of which are outgoing. 



4 Fermion Correlation Functions on Spin-foam 
4.1 Computing Fermion n-point Functions 

Now we consider the correlation functions of a massless Weyl fermion on spin-foam 

(Cii ■ ■ ■ Cj„'?-u„ + i • • • '?i'r. + m)spinfoam 
■= / <i9veT\Af[jf]GRY\^v[jf,ie,gve]GR- f-tdf.i.'X TT '5sL(2.C) eo ) X 

(4.1) 

The Weyl fermion action reads 

-5*6 [6(e), 0(e)' 5'r,6(e)'5'r,/(e)] = [cl(e)5l^b(e)3^./(e)0(e) " ^/(e)ffl,/(e)5reb(e)6(e) 
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where the Dirac operator has the foUowing matrix element 

.t „ , m - _w J 



Pb{e)J(e) = 2Ve grMe)9r,fie) Pf{e),b{e) = 9 (^)gr,b(e) (4.3) 

We employ the standard textbook technique to evaluate the correlation functions. We define a generating functional 

Z{lC,-q\7^)f := V / d5„e TT (^SL(2,C) (Si^eo) IT ^/[-^/iGfl IT ^•"[■?/' *e> S^elcii ' Aj/.Je >< 

= XI / '^Sfe TT <5sL(2,C) (Si'eo ) n ^/ t-^/] Gfl n t-^/ ' *e ' Gi? • A.J/, >< 

X detp[j/,ie,5i'e] exp ^;"j,,[j/,ie,g^e] ?7i,) (4.4) 



Then the correlation function is given by the functional derivative, e.g. the 2-point function (fermion propagator) is 
given by 

(^-i4)spintoam = ^ Z(/C, 7?^ , 7^)/ 

= X / ^5fe TT '^SL(2,C)(5i>eo)n^/[-?/]GflIT^^[-?/'*'='5t'e]Gi? • /7J/,»e >< 

X iPtTi^^sb/'^eiSi^e] det^[j/,ie,5,,e] (4.5) 

A n-point correlation function is given by 

\ / bpmioam 

\ /Spmfoam ^ ^Vv„ + „^r)^Vvi ^Vv„ + „^^^Sl]l^ T,=r,t=0 

= X / '^Sfe TT <5sL(2,C) (S-oeo ) n ^/ t-^/] Gfl n t-^/ ' *e ' Gi? • A.J/, X 

X det p[jf, ie,gve] '^i-'^V Pv^v^+^^.y [jf^^e, gve] ' ' ' Pv^,v„+„^„^ [jf,^e, gve] (4.6) 



where a denotes the permutation (1,2,- •• , n) (ct(1), cr(2), • • • ,cr(n)). The above result shows that the fermion 
correlation function on spin-foams is given by the sum over (disconnected) Feynman diagrams (for free fermion field, 
because here we only consider the interaction between fermions and gravity) 

det p[jf, ie, gve] Y^i-'^Y Pv^.v^^„^^^ [jf^ie, gve] ' ' ' Pv^.v„+,^„y [j/, «e, ff«e] (4.7) 

(T 

on the spin-foam, and the amplitude given by the Feynman diagram depends on the gravitational degree of freedom 
jf,ie,gve, which should also be summed over in the spin- foam model. This result is quite similar to the idea in the 
earlier work [13] about coupling matter Feynman diagrams in spin-foam models. 

The Dirac determinant in Eq.(4.6) has been studied in the previous section. So we are going to find an expression 
for the inverse of the spin- foam Dirac operator -^^^^^^ [j/, iej Sue]- Here we put in a small regulator e > and consider 
the inverse matrix^ (p + e)~^. Since the spin-foam Dirac operator p is anti-Hermitian, the spectrum of operator p + e 



^In case P ^ exists, \img^o{P + ^ = P ^ ■ Even in the ease that P ^ doesn't exist, recaU the Feynman propagator 

^ ^ = 2 ^ (^-^^^ 

If) — m + le p — m + i£ 

We see that the regulator e corresponds to a Feynman regulator for the free quantum field. 
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lies in the right-half complex plane, thus e^*^^^'^-' t e [0, oo) gives a contraction semigroup on a finite-dimensional 
vector space (with the usual norm). Then we have the following strong operator-equation as a consequence of the 
Hille-Yosida theorem [28] : 

/■oo 

dL e-(^+-)^ (4.9) 



which is also known as Schwinger's proper time representation in physics literatures. Given two vertices v,v' 
Jo 

( 1)'^ 
dLe-- E - 



Path,, 



fc! 

(-1)'= 



TdLe-'^ E ^^^'ri(2K.)(5l.5e,„,)(5l,.,5e...)---(5l„,_,5e..')(-l)°P (4.10) 

° Path,,^,,, ' i=l 



where P is nonvanishing only when v,v' are neighboring vertices, X^Path , denotes the sum over all the oriented 
paths from v to v' along the edges of the 2-complex /C, and the number Op denote the number of edges on the path, 
whose orientation are opposite to the path. 

For a Dirac fermion the regulator e should be replaced by iMy_y' + e — iSy^yf^^VyiriQ + e where ^Vi, is the volume 
of a 4-simplex at v and toq is the fermion (bare) mass, in this case the inverse Dirac operator 



/•CO 

' + tM + e)-l, = / dLe-(**^+^)^ 
Jo 



Path„^„/ ' 1=1 

4.2 World-line representation 

There is another representation of {p + £)~\,, in terms of an discretized world-line action®, which is physically inter- 
esting. This representation is obtained by discretizing the exponential e~^^ 



poo 

r°° r 



e " ' 



e ■ 



(4.12) 



In case the number of vertices of /C goes to be large 

L 
n 



dL 



lim 







E [^-""P 



Path" 

















[i - 




L n J 




L n J 



f°° r L 1 r L 

/ dL e~^-^ lim ^ 'J-",^! - 2C«i -ffli«5ei-"i '^i.i.^a - 2^,1,2^2 -ffl^fi 5. 

"J Path" , 



(4.13) 



where Path"_j.„, denotes the set of paths passing through n — 1 vertices except v and v' , y, = 1 if the orientation of 



e coincide with {v,v') and ~ —1 otherwise. We make a change of variable L t-^ -pi to make L has a dimension of 
length / time 



-^0 ^ Path" / > 



(4.14) 



*On Minkowski spacetime, the discussion of the world-line representation of bosonic propagator often can be found in string theory 
textbooks, e.g. [29]. For the fermionic propagator on Minkowski spacetime, the discussion of world-line representation can be found in e.g. 
[30] and the reference therein. 
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where li + I2 + ■ ■ ■ + In = L. Recall that 



9lv9ev' = ni{e)(j' = na{e)a°'{e) 



(4.15) 



where n"(e) is the normalized tangent vector along the edge e at the begin point b{e), then C^^,n"(e) = n"(e) is the 



normalized tangent vector along the edge {v,v') at 6(e). Therefore in the limit n — > 00 



y exp 



JQ 



dl n47(0 W" 7(/) h/ 7(0 



(4.16) 



where na{'^{l)) is the normalized tangent vector along the path ^{l), ¥{^{1)) is the 3-volume of the tetrahedron t{1) 
at 7(0- We then make a change of variable and define t ^ l/L 

-I 



1 

^■p Jo 



e 'p exp 



Jo 



dtin„(7(t))a"(7(t))v-(7(i)) 



(4.17) 



This expression can be written as a gauge-fixed path integral (n"(<) — n"-{j{t)), iT"(t) = a"{j{t))) 
1 



-J^ dL I [De{t)] 



'^'^ \Dj{t)] n 6{eit)-L)e 



ja-/o dt e(t) -jlr/o'dt e(t)n„ (t)<T° (t)^^W 
/Jg p 



(4.18) 



where /^/gV—J' [^7(0] = Spath / nothing but a counting measure. If we define a world-line action (matrix)^ 



2 



dt 



e{t)n^it)a''{t)Vit) + -ee{t) 



(4.20) 



and consider e{t) as a world-line metric. Then the inverse spin-foam Dirac operator is a discretized path integral of 
this world-line action on the spin-foam background 



1 



' + = ^ dL l[De{t)] 



(4.21) 



Obviously, IltGfo 1] ^(^(^) ~ L) is a gauge-fixing for the world-line reparametrization invariance, and the parameter L 
is a world-line Teichmiiller parameter, which is a gauge fixing left-over. 
For Dirac fermion, the fermion world-line action reads 



2 

SD[e,'-f,jf,-ie,9ve] = -^J dt 



e(t)n„(t)7"(<)y(t) + ^ (imo^Vit) + e) e{t) 



(4.22) 



where the discretized version of na{t)j"{t) is Q ^iG^el^Gev' ■ Then the inverse spin-foam Dirac operator is a discretized 
version of world-line path integral 



1 

"'0 



dL / [De{t)] 



(4.23) 



5 Conclusion and Discussion 



We have defined and discussed the fermion quantum field coupled with spin-foam quantum gravity, we have defined 
and explored the properties of the fermion correlation functions on spin-foams, where we have shown there is a spin- 
foam analog of PCT symmetry for spin-foam fermions. The concrete evaluation of the fermion correlations function 



^Actually it can also be considered as a "world-tube" action because of the 3-volume V{t) i.e. 



Sw = -J f e(i) / d^iVdetg 

Jo Jrit) 



na(i)<7°(t) + i£' 



(4.19) 



where we redefine the regulator e' such that e = £'V{t). 
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has also been performed, and the main building bloeks, the Dirac determinant and the inverse Dirac operator, has 
been computed. We have shown that the spin-foam fermion correlation functions can be represented as the Feynman 
diagrams of fermion world-lines imbedded in the spin-foam amplitudes. In this article we have considered only the 
interaction between fermions and gravity, so the Feynman diagram imbedded in the spin-foams are factorized into 
disconnected propagators. We expect that even for interacting matter quantum fields, a similar structure holds, i.e. 
the matter field correlation functions could be represented (at least perturbatively) by Feynman diagrams of the 
interacting fields imbedded in the spin-foam amplitudes. 

In closing, we present a remark about the species doubling problem for lattice fermions. The spin-foam fermions 
are defined with a discrete setting, similar to the fermions in lattice field theory. It is well-known that the formal 
discretization of the fermion action on a lattice suffers the problem of species doubling (see any textbook on lattice 
field theory e.g. [27], see also [31]), while the problem is resolved when the discretized Dirac operator satisfies the 
Ginsparg- Wilson relation. Such a Dirac operator can be constructed from the formal discretization in an overlap 
formulation (overlap fermions) by Neuberger [32], which gives an exact chiral symmetry and anomaly calculation 
(see [31] for a summary). The overlap fermions can also be defined on a curved lattice in the presence of external 
gravitational field [33] , where the chiral symmetry and anomaly calculation are reproduced correctly. Because of these 
results, one might consider if the overlap formulation should be employed to define the spin-foam fermion, instead of 
the formal discretization used in the present work. First of all such an idea could be realized straight-forwardly in 
the formulation of overlap fermion, following the technique for example in [33]. However the overlap formulation (and 
Ginsparg- Wilson relation) make correction for the formal discretized Dirac operator by additional terms proportional 
(and higher order of) the lattice spacing a, which is a semiclassical concept in the context of spin-foam model. Thus 
it seems to us that it is unnatural to implement those corrections fundamentally. But it would be interesting to see if 
those corrections can emerge from some certain semiclassical approximations of spin-foam model. It is not hopeless 
in our opinion for the following reasons: the summing over all the geometries in the spin-foam model make it hopeful 
that the fermion doublers are canceled in a similar way to those on a flat random lattice. Some evidences for this have 
been shown in the context of fermion on Regge gravity [22] , where the fermion propagator are computed numerically 
and display excellent agreement with the continuum field theory. 
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